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ABSTRACT 

The distance to the Orion Nebula Cluster (ONC) is estimated using the rotational 
properties of its low-mass pre main-sequence (PMS) stars. Rotation periods, projected 
equatorial velocities and distance-dependent radius estimates are used to form an 
observational sinz distribution (where i is the axial inclination), which is modelled to 
obtain the distance estimate. A distance of 440 ± 34, pc is found from a sample of 74 
PMS stars with spectral types between G6 and M2, but this falls to 392 ± 32 pc when 
PMS stars with accretion discs are excluded on the basis of their near-infrared excess. 
Since the radii of accreting stars are more uncertain and probably systematically 
underestimated, then this closer distance is preferred. The quoted uncertainties include 
statistical errors and uncertainties due to a number of systematic effects including 
binarity and inclination bias. This method is geometric and independent of stellar 
evolution models, though does rely on the assumption of random axial orientations and 
the Cohen & Kuhi (1979) effective temperature scale for PMS stars. The new distance 
is consistent with, although lower and more precise, than most previous ONC distance 
estimates. A closer ONC distance implies smaller luminosities and an increased age 
based on the positions of PMS stars in the Hcrtzsprung-Russcll diagram. 
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1 INTRODUCTION 

The Orion Nebula cluster (ONC) is among the best studied 
star forming regions. It is perhaps the premier cluster for 
investigating star formation and early stellar evolution be- 
cause it is relatively nearby, very young (< 2Myr) and con- 
tains a large population of stars and brown dwarfs covering 
the entire (sub)stellar mass spectrum (0.01 < M/M© < 30 

- see Hillenbrand 1997; Slesnick, Hillenbrand & Carpenter 
2004). The ONC lies just in front of the dense OMC-1 molec- 
ular cloud, so background contamination of cluster candi- 
dates is small and proper motion studies have successfully 
ascribed membership to almost 1000 stars (Jones & Walker 
1988). The ONC is a focus for understanding the initial stel- 
lar mass function, the evolution of circumstellar discs, early 
stellar angular momentum loss, the influence of high mass 
stars on lower mass siblings and their discs, X-ray activity in 
young stars, the history of star formation and the formation 
of star clusters in general (e.g. see Hillenbrand et al. 1998; 
O'Dell 1998; Lada et al. 2000; Muench et al. 2002; Herbst et 
al. 2002; Preibisch et al. 2005; Huff & Stahler 2006; Shuping 
et al. 2006 among many others). 

The distance to the ONC is quite poorly constrained 

- anywhere between 350 pc and 550 pc, with a most likely 
range of 400-500 pc depending on which techniques are con- 
sidered most reliable (see section 2). The distance is an im- 



portant parameter in determining absolute dimensions, ve- 
locities and mass loss rates in the cluster as well as for esti- 
mating stellar luminosities and hence masses and ages from 
evolutionary tracks. 

In this paper the distance to the ONC is determined 
using the rotational properties of pre main-sequence (PMS) 
stars (axial rotation periods P, and projected equatorial ve- 
locities v sin i, where i is the inclination angle of the spin axis 
to the line of sight), together with estimates of stellar radii. 
From these, the observed distribution of sini can be com- 
pared to a model that assumes random spin axis orientation 
and appropriate contributions from uncertainties in the ob- 
servational parameters. Because the estimated stellar radii 
scale linearly with distance, the distance can be treated as a 
free parameter in the model to optimise the match between 
observed and predicted distributions. 

This technique was developed initially by Hendry, 
O'Dell & Collier Cameron (1993) and used to derive a dis- 
tance of 132 ± 10 pc to the Pleiades, which compares well 
with the currently accepted Pleiades distance from numer- 
ous other techniques (O'Dell, Hendry & Collier Cameron 
1994). Preibisch & Smith (1997) used a variant of the 
method to find a distance to T-Tauri stars in the Taurus star 
forming region. Again, the distance they found of 152±10pc 
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now compares well with distances from Hipparcos parallaxes 
(Bertout, Robichon & Arenou 1999). 

In section 2 previous estimates of the distance to the 
ONC are reviewed. The database of ONC rotation measure- 
ments that were used for this work is discussed in section 3. 
The modelling technique is described in section 4, which also 
discusses systematic effects that must be taken into account. 
The results are presented in section 5 and their robustness 
against variations in the model parameters and assumptions 
is tested. Conclusions and a brief discussion are presented 
in section 6. 



2 PREVIOUS DISTANCE ESTIMATES TO 
THE ONC 

Most previous distance estimates were based on fits to the 
upper main sequence of the ONC in the Hertzsprung-Russell 
(HR) diagram and hampered by colour excesses, variable 
extinction, the nebular background and uncertain binarity. 
Walker (1969) found a distance modulus of 8.37 ± 0.05 
(472 pc) from UBV photometry of relatively unobscured 
stars in the outskirts of the ONC. However, Penston, Hunter 
& O'Neill (1975) obtained optical and near infrared photom- 
etry in a 40 arcminute box around the ONC, determining 
distance moduli of 8.1 ± 0.1 (417pc) from the V,B - V 
colour magnitude diagram (CMD), and only 7.7 (347 pc) 
from the V ,V — I CMD. More recent estimates used the 
large UBVubvy(3 dataset of Warren & Hesser (1978). These 
authors found a distance modulus of 8.42 ±0.24 (483 pc) for 
the Orion OBldl subregion, which is a ~ 30 arcminute di- 
ameter box around the ONC, but excluding the central few 
arcminutes around the central Trapezium cluster of high- 
mass stars. Anthony-Twarog (1982) used the same data with 
somewhat different calibrations to obtain a distance modu- 
lus of 8.19±0.10 (435 pc). Breger, Gehrz & Hackwell (1981) 
used infrared and broad band polarisation measurements to 
exclude objects that showed abnormal reddening or evidence 
for circumstellar material and found that the remaining ob- 
jects have an average distance modulus of 8.0 ±0.1 (400 pc) 
in the V ,B — V CMD. The significant discrepancy between 
for example Walker's (1967) and Breger et al.'s (1981) dis- 
tance estimates points to the fact that there may be sub- 
stantial additional systematic uncertainties present in these 
HR diagram-based results. 

Unfortunately, no light is shed on these varying dis- 
tance estimates by other independent techniques. Only one 
star in the ONC has a direct parallax measurement from 
Hipparcos, yielding a distance of 36llg7 8 pc (Bertout et al. 
1999). A more indirect but frequently cited constraint was 
provided by Genzel et al. (1981). Using an expanding cluster 
parallax method they analysed the proper motions and ra- 
dial velocities of H2O masers in the Kleinmann-Low nebula 
(part of OMC-1), thought to be just behind the ONC, and 
obtained a distance of 480 ± 80 pc. Finally, Stassun et al. 
(2004) found an eclipsing binary just 20 arcminutes south 
of the Trapezium stars. They determined absolute radii and 
effective temperatures for the components and hence esti- 
mated a distance of 419 ± 21 pc from the system photome- 
try. However, Stassun et al. suggest that this binary system 
might be foreground to the ONC because at that distance 
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Figure 1. A Hcrtzsprung- Russell diagram for the 95 stars of the 
ONC rotation sample (see text) assuming an ONC distance of 
470 pc. Model isochrones and evolutionary tracks from Siess et 
al. (2000) are shown. Masses are labelled in solar units. 

the object appears older by a few Myr than the bulk of the 
ONC population. 



3 CONSTRUCTION OF THE DATABASE 

The observed sin i distribution in the ONC was investigated 
using a database founded on the catalogue of ONC axial 
rotation periods constructed by Herbst et al. (2002). This 
catalogue also includes rotation periods measured by Stas- 
sun et al. (1999) and Herbst et al. (2000). From this, all ob- 
jects were selected that had an entry in the Jones & Walker 
(1988) catalogue of proper motions within 15 arcminutes 
of the ONC centre and a spectral type, effective tempera- 
ture (Toff), luminosity (L) and radius (R) estimate listed by 
Hillenbrand (1997). The L and R estimates in Hillenbrand 
(1997) assume an ONC distance of 470 pc. Almost all of 
these stars also have a measurement of their I — K excess 
- i.e. the excess colour over that expected from the photo- 
sphere of a "normal" star of similar spectral type. Measure- 
ments of the equivalent width (EW) of the Can 8542A line, 
which is also diagnostic of accretion, are also available for 
most of the sample from Hillenbrand et al. (1998). 

From this subset, vsini measurements were found in 
the catalogues of Rhode et al. (2001) and Sicilia-Aguilar 
et al. (2005). The latter work was performed at a higher 
resolving power (34 000 versus 21500) and so values from 
Sicilia-Aguilar et al. were preferred for v sin i values less than 
20kms~ 1 . For higher velocities the measurement with the 
best quoted fractional precision was chosen. Objects with 
only upper limits to their v sin i were excluded. Three objects 
(JW 192, JW 381 and JW 710) where Herbst et al. (2002) 
state that the measured periods are highly uncertain were 
also excluded. 

3.1 Measurement uncertainties 

The periods measured for low-mass PMS stars are gener- 
ally very precise and accurate, relying on the co-rotation 
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Table 1. Database of objects considered in this work. The full table is available in electronic form and contains 95 rows. A sample is given 
here to illustrate its content. Column 1 gives the identifying number from Jones &; Walker (1988); columns 2 and 3 list the spectral type 
and T c ff from Hillenbrand (1997); columns 4, 5 and 6 list the vsini, its uncertainty and the source of the rotational velocity membership 
(1 - Rhode et al. [2001], 2 — Sicilia-Aguilar et al. [2005]); column 7 lists the rotational period as given by Herbst et al. (2002); columns 
8, 9 and 10 list the log bolometric luminosity (in solar units), stellar radius (assuming a cluster distance of 470 pc) and A(/ — K) from 
Hillenbrand (1997); column 11 gives the equivalent width of the 8542A Call line from Hillenbrand et al. (1998) [negative indicates an 
emission line]; column 12 lists the derived sini. 



JW SpT T cS vsini Avsini Source P logL/L R/Rq A(I - K) EW[Ca] sini 

(K) (kms" 1 ) (kms" 1 ) (d) (mag) (A) 



3 


K8 


3801.8 


29.8 


4.5 


1 


3.43 


0.24 


3.036 


0.55 


2.0 


0.665 


17 


M4 


3228.4 


19.9 


3.9 


1 


3.14 


-0.20 


2.537 


0.10 


0.8 


0.486 


20 


M3.5 


3296.0 


70.3 


11.6 


1 


0.67 


-0.40 


1.933 


0.10 


0.7 


0.481 


25 


M4.5 


3162.2 


18.0 


4.0 


1 


2.28 


-0.88 


1.209 


0.52 


0.0 


0.671 



of magnetic inhomoheneities (starspots) in the photosphere. 
Providing a sufficient time baseline and sampling frequency 
are obtained, period precisions of better than 1 per cent are 
usually achieved. Studies of objects with periods measured 
at more than one epoch suggest that ~ 90 per cent of period 
measurements have at least this level of precision (Herbst et 
al. 2002) . About 10 per cent of period measurements can be 
catastrophically incorrect. The periods could be too short 
by a factor of two if more than one spotted region on a star 
leads to a "double-humped" light curve that is interpreted 
incorrectly. Alternatively, the measured period can be much 
longer than the true periods because of aliasing with the ~1 
day sampling interval that is present in almost all datasets 
(see discussion in Herbst et al. 2002). Many of these in- 
correct periods have already been weeded out by Herbst et 
al. (2002) and others are suspected (e.g. the three stars ex- 
cluded in the previous subsection), but it is possible that a 
small fraction remain. For now an uncertainty of 1 per cent 
on all the period measurements is assumed. 

The resolution of the rotational velocity studies limits 
the vsini measurements to values of more than llkms -1 
and about 5 kms -1 from Rhode et al. (2001) and Sicilia- 
Aguilar et al. (2005) respectively. Above these limits the 
fractional precision of the v sin i measurements are of order 
10 per cent, though vary from object-to-object. Rhode et al. 
(2001) discuss how problems in calibrating the intrinsic res- 
olution of fibre spectrographs can lead to problematic and 
possibly incorrect values for v sin i when rotational broad- 
ening is close to the detection threshold. Rhode et al. sug- 
gest using v sin i — 13.6 km s _1 as a more plausible detection 
threshold for their data, and that is what is used here. With 
similar reasoning only values of v sin i ^ 10 km s~ are used 
from the work of Sicilia-Aguilar et al. (2005). The sample is 
also filtered of poor quality data with a fractional v sin i un- 
certainty greater than 25 per cent to prevent unnecessarily 
broadening the observed sini distribution. 

Stellar radii are estimated from combining L and T e g 
(Hillenbrand 1997). L comes from extinction-corrected I- 
band magnitudes and bolometric corrections, and assumes a 
distance to the ONC of 470 pc. T c s values are estimated from 
the relationship between spectral type and T ff (with some 
small modifications) proposed by Cohen & Kuhi (1979). 
Hillenbrand (1997) estimates that L is uncertain by about 
0.2 dex due to variability and uncertainties in extinction. 
Random errors in spectral types lead to uncertainties of 
about 0.02 dex in the T cS values. As R oc L" 1/2 T c 2 ff , the 



random errors in R are about 0.11 dex. It is fair to say that 
this uncertainty is itself uncertain, and so values between 
0.07 dex and 0.15 dex are tested in the modelling. Hillen- 
brand (1997) also had concerns that because of difficulties in 
assigning an extinction value, the L and hence R of accreting 
classical T-Tauri stars (CTTS) would be underestimated. As 
the rotational database contains information on the I — K 
excess and EW[Ca], both of which are diagnostic of strong 
accretion, such stars can be optionally excluded. 

The final database contains 95 objects with both ade- 
quate P and vsini measurements and estimates of L, T c g 
and R. These objects are listed in Table [1] (available in elec- 
tronic form only) and the HR diagram is compared with the 
stellar evolutionary models of Siess, Dufour & Forestini in 
Fig. [T] (2000 - the variant with a metallicity of 0.02 and no 
convective overshoot). 

3.2 Selection effects and biases 

The database is subject to a number of selection effects and 
biases (see also the discussion in Herbst & Mundt 2005). 
There is a bias against the inclusion of CTTS and in favour 
of weak-lined T-Tauri stars (WTTS). This is because CTTS 
often show accretion-related, non-periodic variability which 
masks the true rotation period of the star. There is also a 
bias towards objects with shorter rotation periods simply 
because these are easier to measure from datasets with a 
limited time span. Neither of these selection effects will bias 
the observed sin i distribution unless CTTS and more slowly 
spinning stars have non-random spin-axis orientation. 

Of more concern is the observational bias against slow 
rotators and objects with small sin i because of the limited 
sensitivity of the vsini measurements. In principle this can 
be accounted for in our modelling by putting a v sin i cut-off 
in the model. This threshold imposes a smooth roll-off in the 
values of sin i that are capable of detection. The exact shape 
of the roll-off depends on the distribution of true equatorial 
velocities (see section P4~2")) . 

A further consideration is the bias away from low incli- 
nation systems due either to the lack of visibility of starspots 
at equatorial latitudes caused by limb darkening, or to the 
reduced amplitude of spot modulation caused by starspots 
at higher latitudes (see discussion in O'Dell & Hendry 1994). 
The actual latitude distribution of spots on very young stars 
is still debatable. There have been theoretical predictions of 
polar concentrations for spots on rapidly rotating stars with 
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deep convection zones (e.g. Schiissler et al. 1996). However, 
Granzer et al. (2000) predict spots over a wide range of 
latitude for fast-rotating T-Tauri stars, with an equatorial 
concentration in slower rotators. The evidence from Doppler 
imaging of spots on very young T-Tauri stars is mixed. Some 
show spot activity at low latitudes, some at high latitudes 
and other have spots at all latitudes (see Joncour, Bertout & 
Bouvier 1994; Johns-Krull & Hatzes 1997; Neuhauser et al. 
1998). For now I make the assumption that there is some in- 
clination ith) below which starspot modulation is never seen 
(see also O'Dell et al. 1994). This threshold is allowed to 
vary over some plausible range or can be tuned to give the 
best fit to the lower end of the observed sini distribution. 



4 MODELLING THE OBSERVED 
INCLINATION DISTRIBUTION 

4.1 The general approach 

The observational estimates of smi for each star are given 
by 



(sini) c 



(I) 



i (v sin i) 



(l) 



where P b a is the observed period, (v sini) b s is the observed 
projected equatorial velocity, Robs is the estimate of the stel- 
lar radius based on a distance to the ONC of 470 pc, and k is 
a constant which depends on the units used for the various 
quantities. The aim is to model the distribution of (sini) b a 
with a Monte Carlo simulation. To that end I define 

Rttuc ^ 



2tt 



Pi 



rue ^truc 



where the subscript "true" indicates the actual values of 
these three parameters in the absence of measurement un- 
certainties and other systematic effects (see below). I assume 
that P bs is related to Ptruc by 



Pobs =Ptrue(l + S P U). 



(•3) 



where 8p is the fractional uncertainty in the period and U is 
a random number drawn from a Gaussian distribution with 
a mean of zero and unit standard deviation. Similarly, 



(vsini) ob s = «truo sini(l + <5„f7) 



(4) 



where sini is drawn randomly assuming that cosi is dis- 
tributed uniformly between and 1 (i.e. random orientation 
of the spin axes) and 8 V is the fractional uncertainty in the 
v sin i measurements. The true velocity and sin i are split 
into separate factors because the effect of an observational 
lower limit to v sin i can only be modelled properly if a dis- 
tribution of Wtruc is initially specified. For instance a star 
with «true = 50kms _1 cannot be in the Rhode et al. (2001) 
observational sample unless sini > 13.6/50, but clearly this 
threshold changes with ftrue (see section [472]) . 

There is a similar expression relating R h s and iitrue but 
this must also take into account that the true distance, D, 
which is assumed common to all stars in the sample, may 
differ from the 470 pc assumed by Hillenbrand (1997). In 
addition a term must be included that admits the possibility 
that a given star could be an unresolved binary system. In 
which case, L is the sum of two components and Rt TU e could 
be smaller than Robs by a factor, 1 ^ b ^ y/2, for binary 
mass ratios ^ q ^ 1. In this case 



Rob 



/470p< 

Jttruc I ^ — 



6(g) 10 d 



(•») 



where <5i og _R is the uncertainty (in dex) of the logarithmic 
radius estimate (see section 3.1). 

Hence, combining equations 1-5 the model for (sini) b s 
can be expressed as 



(sini)obs = 



D 



470 pc 



6(g) lO'Wfl^s 



(6) 



where Ui, U2, U3 indicate that these are three different ran- 
dom numbers taken from a unit Gaussian distribution. 

To estimate the distance to the cluster I adopt the sim- 
plest, unbiased approach suggested by Hendry et al. (1993). 
If equation [6] is re-expressed as 



(sini)obs = 



D 



470 pc 



(7) 



then taking the average over all the stars in the sample and 
the average over all the trials in a model Monte-Carlo sim- 
ulation, the best estimate for the distance to the cluster is 



D = 470 



((sini) b s ) 
(a) 



pc. 



(8) 



4.2 The true equatorial velocity distribution 

The model distribution for vtmc has some influence on the 
derived cluster distance. The reason is that the lower thresh- 
old for v sin i measurements leads to a higher threshold in 
(sin i) bs for lower values of vtruc - To obtain a good model for 
the distribution of (sinz) b s requires a reasonable description 
of the distribution of ftrue- An incorrect Wtruc distribution 
leads to a systematic error in the distance estimate. 

With the size of sample considered here this error is not 
entirely negligible compared with the statistical uncertain- 
ties in D (see section f5 . 2 [) . Fortunately we can check that 
the assumed Wtruc distribution is reasonable by comparing 
the modelled distribution of v sin i from equation [3] (after 
uncertainties and selection thresholds have been applied) 
with that seen in the data. For instance a flat Utruc distri- 
bution would not explain the observed v sin i distribution in 
the ONC (see section IST2)) . 



4.3 The binary correction factor 

Unresolved binary systems will have systematically overes- 
timated L, overestimated R and hence underestimated sini. 
A neglect of this effect would lead to an underestimated 
distance in equation [S] In older clusters where stars have 
reached the ZAMS, unresolved binaries with q > 0.5 are eas- 
ily detected as they lie significantly above the cluster single 
stars locus in the HR diagram. This becomes impractical 
in young clusters like the ONC where there may be an age 
spread that is a significant fraction of the mean age, and 
where differential extinction and variability cause significant 
spreads in any CMD. 

I have calculated a correction to R that should be ap- 
plied to an unresolved binary system of total luminosity L 
and mass ratio q. The calculation uses a model of L versus 
T e ff calculated at a range of masses. For this paper I have 
used the models of Siess et al. (2000) and Baraffe et al. (2002 
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Figure 2. The multiplicative binary correction factor as a func- 
tion of mass ratio for ages of 1 Myr or 3 Myr, for —0.5 < 
log(L / Lq) < 1.0 and for two different evolutionary models (Siess 
et al. 2000 - SOO and Baraffe et al. 2002 - B02). 



- the variant with a mixing length equal to the pressure scale 
height). For an assumed age and L, an initial guess (an up- 
per limit) for mi is taken straight from the isochrone and 
an initial guess for is qm\ . The value of mi is iteratively 
reduced until the total luminosity of the two components 
equals L. At this point the actual T a g of the primary and 
secondary can be calculated as well as a flux-averaged T e g 
that is assumed to correspond to the observed spectral type. 
The multiplicative overestimate of the radius is given by 



(9) 



In Fig. [2] I show AR(q) for -0.5 < logL/L© < 1 and 
for ages of 1 or 3 Myr that cover the majority of stars in the 
ONC rotation sample. At low luminosities the Baraffe et 
al. (2002) models are used because the Siess et al. (2000) 
models do not extend to low masses (M < 0.1 Mq), on 
the other hand the Baraffe et al. models do not extend to 
high luminosities as they are limited to masses less than 
1.4 Mq. At luminosities where the models overlap there is 
good agreement. These curves demonstrate that AR(q) is 
quite insensitive to changes in age, but does behave differ- 
ently for older stars with high luminosity. This is caused by 
a sharply changing slope in the luminosity-mass relationship 
as stars develop a radiative core. This occurs at about 10 Lq 
at 3 Myr, but at 30 Lq for ages of ^ 1 Myr. This boundary 
more-or-less coincides with the envelope of the sample stars 
considered here (see Fig. [1} and so a single AR(q) curve was 
adopted that is representative for most of the sample (the 
Siess et al. model with logL/L© = 0.5 at 3 Myr), which is 
given the label "rl" in sect ion [5T3l To check the sensitivity of 
the results to this assumption the most deviant AR(q) curve 
in Fig. [5] (the Siess et al. model with logL/L0 = 1.0 and at 
3 Myr) is also tested in section [5T31 and given the label "r2". 
The difference in derived distance between using models rl 
and r2 turns out to be small (2 per cent), but the neglect of 
binaries altogether would lead to a distance underestimate 
of ~ 10 per cent. 

The relationship between AR(q) and b(q) in equation[6] 



is set by assuming that a fraction / of star-systems are bi- 
naries, with q drawn randomly from a defined distribution. 
The appropriate value of / and the q distribution can be 
inferred from observational studies of the ONC. 

Using different high spatial resolution surveys Prosser 
et al. (1994), Padgett, Strom & Chez (1997), Petr et al. 
(1998), Simon, Close & Beck (1999) and Kohler et al. (2006) 
have concluded that the binary frequency in the range 100- 
1000 au is similar to that of field stars and probably declines 
towards lower masses in a similar way. Details for closer bi- 
nary systems are yet to be published, but preliminary results 
suggest that the close binary frequency is also similar to that 
of field stars (Stassun & Mathieu 2006). 

I make the assumption that binary properties in the 
ONC are similar to that of the field. Duquennoy & Mayor 
(1991) show that the overall binary frequency for solar-mass 
stars is 57 per cent and at distances of ~ 450 pc, the vast 
majority of these would be unresolved binaries in the ONC 
sample. The binary frequency declines towards 30-40 per 
cent in M dwarfs with 0.3-0.5 Mq (Fischer & Marcy 1992). 
The mass ratio distribution appears reasonably uniform, at 
least over the range 0.5 < q ^ 1 that is most important for 
the phenomena considered here. The stars in our rotational 
sample cover a mass range of 0.2 < M/Mq < 2.5 according 
to the PMS tracks of Siess et al. (2000; see Fig. [l}. Rather 
than complicate things further, the binary frequency in the 
ONC is assumed to be 50 per cent, mass-independendent 
and q is given a uniform distribution. The sensitivity of the 
results to these assumptions is tested by allowing the binary 
frequency to vary by ±20 per cent in the simulations (see 
section I5.3|l . 

It is worth stressing that one could just use an aver- 
age value of 6(g) as a multiplicative correction factor to the 
distance estimate in equation [5] However, I calculate indi- 
vidual b values for each star in the Monte Carlo simulations 
because binarity introduces an extra spread in the model a 
distribution which should be taken into account when decid- 
ing whether the model is a reasonable fit to the data and, 
for instance, in deciding on an appropriate value for ith (see 
section 3.2). 

As a final remark on binary systems, I note that the 
vsini measurements of very short period (less than a few 
days) binaries could be artificially increased by orbital mo- 
tion because of the lengthy exposure times used to obtain the 
spectra (up to a few hours for the Rhode et al. [2001] mea- 
surements). This would increase (sini) b s for the affected 
objects, leading to a distance overestimate. I have neglected 
this possibility in the modelling, but as only a very small 
minority of binary systems are likely to have such short pe- 
riods and be observed at unfavourable phases, this seems 
reasonable. 



4.4 Recap of assumptions and the step-by-step 
method 

Having described the general methods I use to calculate a 
model (sini) bs distribution, the following summarises the 
specific assumptions and steps in the modelling process. 

(i) Choose the sample. Stars are selected from the ro- 
tation sample on the basis of having a minimum observed 
projected equatorial velocity («sini) b s (13.6 km s^ 1 and 
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lOkms -1 for the Rhode et al. [2001] and Sicilia-Aguilar 
[2005] observations respectively - see section 3.1) and can 
be further restricted to a range of temperature and whether 
they show signs of accretion or not (see section [5~H . 

(ii) A model is chosen for the true equatorial velocity 
(utrue) distribution. Then, assuming that the rotation axes 
are randomly orientated, a set (typically 10 4 trials per object 
in the sample) of model usini values are randomly gener- 
ated and perturbed according to equation 2] using the frac- 
tional uncertainties in (v sin i) bs as the S v values. The model 
^truc distribution is adjusted until a satisfactory match to 
the (usini) bs distribution is indicated by a Kolmogorov- 
Smirnov (K-S) test between the observed and modelled cu- 
mulative distribution functions (CDFs). 

(iii) The trial values of sin i are used to calculate the a 
values in equation [7] using S v , fractional errors in the period 
5p = 0.01, and assuming a suitable value for the fractional 
error in the log stellar radius (<5iogfl). A fraction / of the 
trials are randomly assigned a binary status and the appro- 
priate correction factor to the observed radius, b(q), is ap- 
plied assuming that the binary mass ratio (q) is uniformly 
distributed. 

(iv) The results of the trials pass through a filter which 
only allows those trials to proceed which have v sin i greater 
than the observational thresholds and axial inclinations 
greater than a user-defined threshold, i t h- 

(v) At this stage the average value of (sini) b s is divided 
by the average value of a in order to find D (equation [8]) . 
At the same time a K-S test is made between the CDFs of 
a(D/470) and (sinz) b s . This is used to decide whether the 
model is a reasonable description of the data (the sample is 
too small to consider x 2 tests). 

(vi) Statistical uncertainties in the results are estimated 
by generating many fake datasets of the appropriate size. 
The (sini) bs values are drawn randomly using equation[7]at 
the estimated distance and these are modelled in the same 
way as the data. The standard deviation of the distance 
estimates is used as a statistical uncertainty. 

(vii) The sensitivity of the results to the model assump- 
tions and parameters are tested using different subsamples 
of data (e.g. with or without signs of accretion), different 
values of i t h, /, 8\ og R or different models for AR(q) and the 
^truc distribution (see section I5~31 . 



5 RESULTS 

5.1 The observed sini values 

In Fig. |3] (sin i) bs versus T c ff is plotted for the sample. The 
data have been distinguished on the basis of whether their 
near infrared excesses, A (7 — K), are larger or smaller than 
0.3. This division approximately separates stars with and 
without active accretion and thus the two subsets are la- 
belled as CTTS and WTTS respectively. Hillenbrand et al. 
(1998) show that because of physical effects (disc inclina- 
tion, inner holes etc) and also observational uncertainties, 
this threshold is fuzzy and not capable of establishing or 
excluding accretion in all cases. Unfortunately, mid-IR mea- 
surements which are much more sensitive are only available 
for a small fraction of the sample (see Rebull et al. 2006). 
An alternative is to use EW[Ca], which goes into emission 
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Figure 3. The observed value of sini (from equation 1) versus the 
effective temperature. CTTS/WTTS are indicated by different 
symbols. The lines show low order spline fits as a function of T e g 
to (a) all the data, (b) the CTTS only (A(I - K) > 0.3), (c) the 
WTTS (A(J - K) < 0.3). 
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Figure 4. Similar to Fig. |3] but here the WTTS and CTTS are 
classified according to whether the equivalent width of the Can 
8542A line is more or less than lA (see Hillenbrand et al. 1998). 

for strongly accreting objects. WTTS are generally found to 
have EW[Ca]> lA and Fig.|4]repeats Fig.[3]using this crite- 
ria to separate WTTS and CTTS. Although crude, these di- 
agnostics should be sufficient to identify cases where broad- 
band colours, magnitudes and hence estimated luminosities 
and radii are likely to be affected, which is the issue here. 

Two sample selection issues are highlighted by Figs. [3] 
and [4] The first is that stars with a near infrared excess or 
EW[Ca]< lA have a larger (sini) b s on average than those 
without. To illustrate this, smooth spline fits as a function 
of T off are plotted for all the data and for the CTTS/WTTS 
subsets separately. The effect is small at cool temperatures 
but increases to ~ 30 per cent in hotter stars. For stars with 
T cf f > 3499 K (see below) the average sini is 0.78 ± 0.04, 
but the averages for the WTTS samples are 0.69 ± 0.04 for 
34 stars with A(J — K) < 0.3 and 0.67 ± 0.04 for 44 stars 
with EW[Ca]> lA. As the estimated distance is linearly de- 
pendent on (sini) ba, this bias is of concern. There are no 
physical reasons to suppose that CTTS have larger inclina- 
tion angles, but there are good reasons to suppose that the 
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radii of CTTS have been underestimated (leading to a sin i 
overestimate), either because of a systematic underestima- 
tion of the extinction and luminosity (see Hillenbrand 1997) 
or because of difficulties in determining their temperatures 
(see below and section 6). 

The second issue is that the average (sini) b s appears 
to get smaller for T c g < 3499 K (corresponding to spectral 
types cooler than M2) in both WTTS and CTTS. The aver- 
age (sini) obs for WTTS (A(i-it) < 0.3) with T cff < 3499 K 
falls to only 0.53 ± 0.04. There are a number of possible 
causes for this, (i) Perhaps the binary fraction is a strong 
function of mass (and therefore T c g). A larger binary frac- 
tion, or tendency towards high q systems will increase the 
average i? bs and hence reduce the average (sini) b s - The 
evidence from field stars is that the binary fraction falls sig- 
nificantly towards lower masses, so this explanation seems 
unlikely, (ii) The observed radius is derived from L and a 
T e flf which is derived from the Cohen & Kuhi (1979) main- 
sequence relationship between T e s and spectral type. For 
warmer stars (T e a ^ 3500 K) this relationship is relatively 
uncontroversial in the sense that alternative scales (for in- 
stance Kenyon & Hartmann 1995; Leggett et al. 1996) differ 
by less than 2 per cent and not in a systematic way. For 
cooler stars the relationship is more uncertain and espe- 
cially so for very cool PMS stars which have significantly 
lower gravities than main-sequence stars of similar spectral 
type. The luminosity estimates are also reliant on an ac- 
curate assessment of the extinction which requires intrin- 
sic broadband colours for very cool PMS stars. Hillenbrand 
(1997) used the intrinsic colours of main-sequence stars tab- 
ulated by Bessell & Brett (1988), but very cool PMS stars 
may have more giant-like colours. Dwarfs are bluer than gi- 
ants by ~ 0.3 mag in V — I for spectral types M3-M5, which 
would lead to a larger extinction estimate by ~ 0.5 mag 
(at 7). This gives a larger luminosity estimate by 0.2 dex, 
a larger radius estimate by 0.1 dex and hence a decreased 
sini estimate by a factor of 1.26 that would almost entirely 
explain the trend seen. 

For the reasons above the analysis here is initially re- 
stricted to the sini distribution of those 34 objects in Ta- 
bled] with A(J - K) < 0.3 and T eS > 3499 K. This sample 
is labelled WTTS1 in what follows. Two alternate samples 
are also considered, namely 44 objects with T ff > 3499 K 
and EW[Ca]> lA (labelled WTTS2) and all 74 objects with 
T c ft > 3499 K, regardless of their accretion status (labelled 
ALL). 

5.2 The vtrue distribution 

For the reasons explained in section 14.21 a reasonable model 
of the vtruc distribution is required and can be constrained 
by comparing the simulated v sin i distribution (from many 
trials) with that observed. 

In Figure two of these comparisons are shown. A 
formal statistical comparison has been made using a 1- 
dimensional K-S test of the CDFs. The models shown have 
ith = 30° , but the model v sin i distribution does not depend 
strongly on this parameter. 

After experimenting with a variety of simple analyt- 
ical forms, I find that a two component v^uc model (la- 
belled vl in Fig. [5} provides a good (though not neces- 
sarily unique) fit to the data, with 80 per cent of stars 



following an exponentially decaying distribution between 
10 < utruc < 120 kms" 1 , with a decay constant of 20 kms -1 , 
and the remaining 20 per cent being uniformly distributed 
between 10 < utruc < 120kms _1 . The K-S test indicates a 
probability of 85 per cent for the null hypothesis that the 
the data is drawn from the model distribution. On the other 
hand, the K-S test shows that a flat Utruc distribution (la- 
belled v2 in Fig. [SJ can be excluded with 99.999 per cent 
confidence. 

5.3 Distance estimates 

Following the procedure in section 4.4, the distributions of a 
and (sini) bs were calculated and used to estimate the clus- 
ter distance from equation 8. A K-S test comparison was 
performed between the CDFs of a(D/470) and (sini) b s to 
test whether the model is a reasonable representation of the 
data. Statistical uncertainties were estimated by generat- 
ing fake datasets of the same size, with Utruo and binary 
properties drawn from the distributions used in each model. 
Observational errors for the fake datasets matched those in 
the real dataset considered here. The fake datasets passed 
through the same analysis procedure and the standard de- 
viations of the distance estimates (from 300 fake datasets) 
were used as estimates of the statistical error. 

The results are collected in Table [2] and Fig. [6] demon- 
strates the sensitivity of the model distribution to changes 
in the various parameters in the form of both the differ- 
ential and cumulative distributions of a(D/470) compared 
with those of the data. 

Analysis number 1 from Table [5] is used as a baseline 
for the discussion as this uses the parameters that could be 
considered to most likely represent the modelled ONC pop- 
ulation. The average value of (sini) b s is significantly lower 
than would be expected for randomly inclined axes in agree- 
ment with the analysis of Rhode et al. (2001) based on a 
smaller sample. The simplest explanation is that the ONC is 
much closer than the 470 pc assumed by Hillenbrand (1997). 
Indeed D = 470 pc is rejected with 97 per cent confidence. I 
find a distance of 392 pc gives a much better representation 
of the data as may be judged from Figs. |6^, and [BJj. The 
statistical uncertainty is ±23 pc which compares well with 
simulations and uncertainty estimates presented by O'Dell 
et al. (1993) for similar sized samples. 

In addition to the statistical uncertainties I have tried 
other models to test the sensitivity of the derived distance 
to variations in model parameters that are only partially 
constrained. Fig. [5] shows examples of the model distribu- 
tions at a fixed distance of 392 pc in order to highlight these 
differences. The reader should note that in all cases the vari- 
ations in the intrinsic a distribution can be compensated for 
by changes in the cluster distance to yield a reasonable value 
for P(K-S) (these are the results quoted in Table [2}. There- 
fore the observed (sini) (, s distribution is currently incapable 
of constraining some of the more interesting parameters like 
ith, due to the small sample size and experimental errors. 

Figs. and[B}i show how changes of 15° < i t h < 45° 
affect the modelled distribution at a fixed distance of 392 pc. 
The changes are small and easily compensated for by mod- 
est changes in the cluster distance of only ±4.0 per cent. 
This range seems a reasonable estimate of the uncertainties 
engendered by the inclination bias - if ith < 15° this would 
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Figure 5. A comparison of the observed usini distribution with Monte-Carlo simulations based on two different models of the vtmc 
distributions, vl and v2 (described in the text). The left panel shows the (binned) distributions and the right panel shows the normalised 
cumulative distributions used to perform a Kolmogorov-Smirnov test. The discontinuity in the model distributions at usini = 13.6 km s — 1 
is due to the modelled sample being divided between objects with detectable usini thresholds at 13.6kms — 1 and lOkms - . 



Table 2. The results of the distance analysis. The columns are: (1) An analysis identification number; (2) the analysed sample and 
number of objects in the sample - WTTS1 indicates objects A(J - K) < 0.3, WTTS2 indicates objects EW[Ca]> lA (see section 1570 . 
ALL indicates no A(I — K) restriction; (3) identifies which AR(q) model was used, rl or r2 (see section 14.31 1; (4) identifies which 
vtruc distribution was assumed, vl or v2 (see section [5.2H ; (5) the assumed value of i t h (see section 3.2); (6) the assumed value of the 
binary fraction / (see section [4.31 1; (7) the assumed value of uncertainty in the observed (log) radius (see section 3.1); (8) the average 
value of (sini) t, s ; (9) the average value of a; (10) the distance implied by equation 8 and its 1-sigma statistical uncertainty; (11) the 
Kolmogorov-Smirnov null-hypothesis probability (the probability that the observed distribution is not drawn from the model distribution 
is 1 - P(K-S)). 



Analysis No. 


Sample (n) 


AR(q) 


Wtrue 


kh 


bfrac 


<5io g fl 


((sini) obs ) 


(a) 


D 


P(K-S) 






model 


model 


degrees 


/ 


(dex) 






(PC) 




1 


WTTS1 (34) 


rl 


vl 


30 


0.5 


0.11 


0.691 


0.828 


392 ± 23 


0.66 


2 


WTTS1 (34) 


rl 


vl 


15 


0.5 


0.11 


0.691 


0.802 


405 ± 24 


0.75 


3 


WTTS1 (34) 


rl 


vl 


45 


0.5 


0.11 


0.691 


0.870 


373 ± 22 


0.61 


4 


WTTS1 (34) 


rl 


vl 


30 


0.3 


0.11 


0.691 


0.858 


378 ± 22 


0.67 


5 


WTTS1 (34) 


rl 


vl 


30 


0.7 


0.11 


0.691 


0.798 


407 ± 24 


0.63 


6 


WTTS1 (34) 


rl 


vl 


30 


0.5 


0.07 


0.691 


0.812 


400 ± 18 


0.64 


7 


WTTS1 (34) 


rl 


vl 


30 


0.5 


0.15 


0.691 


0.851 


382 ± 28 


0.41 


8 


WTTS1 (34) 


rl 


v2 


30 


0.5 


0.11 


0.691 


0.815 


398 ± 24 


0.64 


9 


WTTS1 (34) 


r2 


vl 


30 


0.5 


0.11 


0.691 


0.849 


383 ± 22 


0.68 


10 


WTTS2 (44) 


rl 


vl 


30 


0.5 


0.11 


0.671 


0.829 


381 ± 20 


0.66 


11 


ALL (74) 


rl 


vl 


30 


0.5 


0.11 


0.780 


0.834 


440 ± 19 


0.48 



make little further difference to the distance estimate and 
if i t h > 45° then there would be difficulties in measuring 
rotation periods in a large fraction of WTTS. 

Figs. [6J3 and [6J show the effects of changing the bi- 
nary frequency between 0.3 < / < 0.7 at a fixed distance. 
Again the small differences can be compensated for by small 
changes (±3.7 per cent) in distance - larger / leads to larger 
distances. 

Figs. |6g and[6h show that changing the assumed value 
of Sio g R in the plausible range 0.07 < <5i og fl < 0.15 sig- 
nificantly changes the shape of the distribution at a fixed 
distance. This can be partially compensated for by chang- 
ing the distance by ±2.3 per cent. The broader distribution 
is a poorer fit to the cumulative distribution function but 
cannot be rejected. 

Figs.|U and[S) show that even taking implausibly large 
variations in the v tIU c and AR(q) distributions (see sec- 
tions [5]2] and |4]3] respectively) results in very small changes 



in the model that can be compensated for by changes in the 
cluster distance of only ±1,9 per cent. 

For completeness, I show the effects of considering 
the sample of 74 objects (ALL) including both CTTS and 
WTTS (but all with spectral types of M2 or hotter) and also 
a sample of WTTS defined on the basis of EW[Ca] (WTTS2 
- see section[S]TJ . As expected from Figs.|3]and|3]and the dis- 
cussion in section [5Tl the distribution for the entire sample is 
skewed to higher (sini) b s values, resulting in a significantly 
larger distance estimate with slightly better statistical pre- 
cision (440 ± 19 pc). The sample of WTTS based on EW[Ca] 
is larger than sample WTTS1 but filters out a few more of 
the objects with large (sini) b s values. Hence the deduced 
distance based on sample WTTS2 is marginally closer than 
that based on sample WTTS1. 

The "best" final ONC distance estimate is 392 pc. The 
total uncertainty is 32 pc, estimated from quadratic sum of 
the statistical uncertainty and the distance uncertainties due 
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Figure 6. The observed (sini) t, s distribution compared with various model distributions of the form a(D/470 pc). The left hand panels 
show the differential distributions while the right hand panels show the normalised cumulative distributions, (a, b) Shows the effects of 
distance on the data-model comparison. The three models shown correspons to model 1 in Table |51 at D = 392 pc, model 1 at D = 415 pc 
(corresponding to a 1-sigma error) and model 1 at D = 470 pc. (c, d) Shows the effects of varying i t h by displaying models 1, 2 and 
3 (*th = 30°, 15° and 45°) at a common D = 392 pc. (e,f) Shows the effects of varying the binary fraction / by displaying models 
1, 4 and 5 (/ = 0.5, 0.3 and 0.7) at a common D = 392 pc. (g, h) Shows the effects of varing 5i og i? by displaying models 1, 6 and 7 
(^logfl = 0.11 dex, 0.07 dex and 0.15dex) at a common D = 392 pc. 



© 2004 RAS, MNRAS 000.ITHT21 



10 R.D. Jeffries 



12 
10 

8 
6 
4 
2 


20 
15 
10 




v2, rl,D=392 pc 
vl,rl,D=392 pc 
vl,r2, D=392 pc 



( sin i) obs 



(sin i) 



□ 



WTTS2, D=381 pc 
ALL, D=440 pc 





1 

0.8 
0.6 
0.4 
0.2 



WTTS2, D=381 pc 
ALL, D=440 pc 




( sin i)ni 



Figure 5. continued: (i, j) Shows the effects of choosing different v truc (vl, v2) and AR(q) (rl, r2) distributions by displaying models 
1, 8 and 9 at a common D = 392 pc. (k, 1) Shows a model which best fits a sample of WTTS defined using EW[Ca] (see section |5.1| | 
(model 10 at D = 381 pc) and a model which best fits a sample containing both CTTS and WTTS (model 11 at D = 440 pc). 



to the parameter variations discussed above. This of course 
assumes these uncertainties are independent, ignores any 
systematic problem with the spectral type-effective temper- 
ature relationship, and assumes that it is correct to exclude 
CTTS from the analysis. The equivalent figure including the 
CTTS would be 440 ± 34 pc, where the uncertainty in this 
case is dominated by systematic effects. 



6 DISCUSSION AND CONCLUSIONS 

The observed sini distribution of a group of WTTS (and 
CTTS) in the ONC has been modelled. The sini values are 
linearly dependent on the cluster distance and from this a 
distance to the ONC of 392 ± 32 pc is derived or 440 ± 34 pc 
if CTTS are included in the sample. The uncertainties in- 
clude 1-sigma statistical uncertainties and the contributions 
from uncertainties in a number of other partially constrained 
model parameters such as the binary frequency, bias against 
the inclusion of low inclination objects, the form of the true 
equatorial velocity distribution and uncertanties in the stel- 
lar radius estimates. With the present analysis and observa- 
tional sample (34 WTTS), statistical and systematic uncer- 
tainties are of equal importance. 

The observed sini values depend linearly on the esti- 
mated radii, which are difficult to determine for PMS stars. 
Hillenbrand (1997) discusses these difficulties and concluded 
that extinctions, luminosities and hence radii may be under- 
estimated for accreting CTTS, resulting in overestimates of 
sin i and hence overestimates of the distance using the anal- 
ysis presented in this paper. This effect is readily apparent 
in the data, even when using imperfect discriminators of ac- 



cretion activity such as the I — K excess and EW of the 
Can 8542A line. A more careful analysis could use mid-IR 
measurements to exclude (accretion) discs, but for now it 
seems wisest to discard the suspected CTTS and adopt the 
closer of the two distance estimates above. 

A further issue is the validity of the spectral- 
type/effective temperature relationship. As stated in sec- 
tion [STT1 the relationship of Cohen & Kuhi (1979) is uncon- 
troversial in that it is consistent with several later studies. 
However, these are also fundamentally based on the atmo- 
spheres, angular diameters and bolometric luminosities of 
cool main sequence stars. Young PMS stars are likely to have 
different atmospheres and the fact that rotational modula- 
tion has been measured in our sample stars means that they 
must be covered by cool spotted regions. Assigning a single 
temperature to such magnetically active stars may be diffi- 
cult. These difficulties could be worse in CTTS with their 
attendant discs, veiling and accretion hotspots and this is a 
further reason for excluding them from the analysis. As yet, 
there is no compelling evidence for a change in the spectral- 
type/temperature relationship for non-accreting PMS stars, 
but as deduced distances will scale as T^ s this issue should 
be borne in mind. 

The deduced ONC distance is smaller than most previ- 
ous results, although generally agrees within the error bars. 
If proper account of systematic uncertainties is included, the 
distance derived here is also probably more accurate than 
most previous determinations. It is worth noting however 
the excellent agreement with the 419 ± 21 pc determined by 
Stassun et al (2004) from an eclipsing binary that they con- 
sidered to be a possible ONC member. 

To recover the more usually used longer ONC distance 
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of 450 pc or 470 pc (e.g. Hillenbrand 1997; Preibisch et al. 
2005) from the rotational data would require a model with 
a very high binary fraction (/ ~ 1) combined with a sev- 
eral per cent change in the relationship between spectral 
type and temperature for warm K-M2 ONC stars. A closer 
distance reduces deduced luminosities in the ONC by a cor- 
responding factor. Perhaps more importantly, ages deduced 
from the positions of PMS stars in the HR diagram are in- 
creased by a factor of ~ 1.7, because L oc t~ 2 ^ 3 on fully 
convective PMS tracks. 

The distance determined here is independent of stellar 
evolutionary models - an important source of systematic er- 
ror in HR diagram-based estimates. However, it does hinge 
on the assumption that the rotation axes of PMS stars are 
randomly oriented in space. The observed sini distribution 
is certainly consistent with this hypothesis, but is also suffi- 
ciently broadened by experimental uncertainties that other 
hypotheses cannot be rejected. For example, taking the pa- 
rameters of model 1 in Tableland using a fixed inclination 
i — 60° results in a distance of 396 pc and P(K-S)= 0.48. 

The theoretical and observational justification for ran- 
dom axial orientation is mixed. In the quasi-static picture 
of star formation, cloud collapse occurs preferentially paral- 
lel to magnetic field lines. Some polarization studies suggest 
that disc (and therefore presumably rotation) axes could be 
aligned with the large scale magnetic fields in star forming 
regions (e.g. Vink et al. 2005). On the other hand, imag- 
ing of CTTS discs in the Taurus-Auriga region reveals a 
random orientation with respect to the large scale magnetic 
field, in accord with a more dynamical star formation picture 
(Menard & Duchene 2004). 

As a final remark, previous applications of the rota- 
tionally based technique to the Pleiades cluster and Taurus- 
Auriga association have yielded distances which turn out to 
be in good agreement with later precise parallax-based mea- 
surements (see O'Dell et al. 1994; Preibisch & Smith 1997). 
This tends to argue that the underlying assumptions of the 
method (random axial orientations and a dwarf-like temper- 
ature scale for non-accreting PMS stars) are reasonable. A 
subsequent precise trigonometric parallax for the ONC (e.g. 
using the Hubble Space Telescope) could be used to test the 
robustness of these model assumptions in more detail. 
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